Solar Physics 
DOI: 10.1007/' 



Magnetic Field Configuration Models and 
Reconstruction Methods for Interplanetary Coronal 
Mass Ejections 

N. Al-Haddad 1 2 , T. Nieves-Chinchilla 3 4 , 

N. P. Savani 5 6 , C. Mostl 7 8 9 , 

K. Marubashi 10 11 , M. A. Hidalgo 12 , 

I. I. Roussev 2 , S. Poedts 1 , C. J.Farrugia 13 



© Springer •••• 

Abstract This study aims to provide a reference to different magnetic field 
models and reconstruction methods for interplanetary coronal mass ejections 
(ICMEs). In order to understand the differences in the outputs of those models 
and codes, we analyze 59 events from the Coordinated Data Analysis Workshop 
(CD AW) list, using four different magnetic field models and reconstruction tech- 
niques; force-free fitting (Goldstein, 1983; Burlaga, 1988; Lepping, Burlaga, and 
Jones, f990), magnetostatic reconstruction using a numerical solution to the 
Grad-Shafranov equation (Hu and Sonnerup, 2001), fitting to a self-similarly 
expanding cylindrical configuration (Marubashi and Lepping, 2007) and ellip- 
tical, non- force free fitting (Hidalgo, 2003). The resulting parameters of the 
reconstructions for the 59 events are compared statistically, as well as in selected 
case studies. The ability of a method to fit or reconstruct an event is found to 
vary greatly: the Grad-Shafranov reconstruction is successful for most magnetic 
clouds (MCs) but for less than 10% of the non-MC ICMEs; the other three 



1 Centrum voor Plasma- Astrofysica, Katholieke Universiteit 
Leuven, Celestijncnlaan 200B, 3001 Leuven, Belgium 
2 Institute for Astronomy, University of Hawaii, 2680 
Woodlawn Dr., Honolulu, HI 96822, USA. 
3 Heliospheric Physics Lab. GSFC-NASA, Greenbelt, MD. 
USA. 4 IACS-CUA, Washington, DC, USA. 5 UCAR, 
Boulder, CO, USA 6 NASA Goddard Space Flight Center, 
Greenbelt, MD, USA. 7 Space Science Laboratory, University 
of California, Berkeley, CA, USA. 

8 Space Research Institute, Austrian Academy of Sciences, 
Graz 8042, Austria. 

9 Kanzclh6hc Observatory-IGAM, Institute of Physics, 
University of Graz, Universitatsplatz 5, A-8010, Graz, 
Austria. 

1(, Korca Astronomy and Space Science Institute, Dacjcon, 
305-348, Republic of Korea. n 660-73 Mizuno, Sayama, 
Saitama 3350-1317, Japan. 

12 SRG-UAH, Alcala de Henares, Madrid, Spain. 
13 Space Science Center and Department of Physics, 
University of New Hampshire, Durham, NH, USA. 



SOLA: Comparison_PaperAfter_revs.tex; 1 October 2012; 0:39; p. 1 



Al-Haddad et al 



methods provide a successful fit for more than 65% of all events. The differences 
between the reconstruction and fitting methods are discussed, and suggestions 
are proposed as to how to reduce them. We find that the magnitude of the axial 
field is relatively consistent across models but not the orientation of the axis 
of the ejecta. We also find that there are a few cases for which different signs 
of the magnetic helicity are found for the same event when we do not fix the 
boundaries, illustrating that this simplest of parameters is not necessarily always 
well constrained by fitting and reconstruction models. Finally, we look at three 
unique cases in depth to provide a comprehensive idea of the different aspects 
of how the fitting and reconstruction codes work. 

Keywords: Sun: corona — Sun: coronal mass ejections (CMEs) 
1. Introduction 

Coronal mass ejections (CMEs) are large-scale eruptions of plasma and magnetic 
flux with a typical size of 0.25 AU at 1 AU (e.g., see Bothmer and Schwenn, 1998). 
In the interplanetary space, they typically move at supersonic speeds between 
200 and 2000 km s _1 , propagating from the solar corona to 1 AU in a few days. 
They can interact with Earth's magnetosphere resulting in geomagnetic storms 
(Gold, 1959). Once CMEs drive through and interact with the interplanetary 
medium, they are often referred to as interplanetary CMEs (ICMEs). ICMEs, 
measured in situ, for example at LI, may be composed of a fast forward shock, 
a dense sheath and ejecta material. Slow CMEs typically do not drive a shock. 
Hereafter, we use the term ICME once it is propagating in the interplanetary 
medium after leaving the corona (more or less after 10-20 Rq). 

ICMEs, measured in situ, are typically characterized by low proton tempera- 
ture, the presence of bi-directionally streaming electrons, unusual charge states of 
oxygen and iron and various magnetic signatures (e.g. Lynch et al., Zurbuchen 
and Richardson, Richardson and Cane, 2003, 2006, 2010). However, there is 
no single characteristic that is consistently observed. Lists of ICMEs measured 
by ACE and STEREO are maintained by Jian et al. (2006), Lepping et al. 
(2006), Richardson and Cane (2010) and Gopalswamy et al. (2008), among 
others. ICMEs have been classified into three different classes based on their 
magnetic field and plasma properties (Zurbuchen and Richardson, 2006): 

1. Magnetic clouds (MCs): characterized by low plasma beta, low proton tem- 
perature and featuring a strong magnetic field with a smooth rotation (Burlaga 
et al, 1981). In order to explain the last characteristic, Goldstein (1983) sug- 
gested that the magnetic field in a MC can be described by a force-free flux 
rope structure, where V x B = aB. Burlaga (1988) took a constant a, in which 
case the flux rope satisfies the Lundquist solution (Lundquist, 1950) and can be 
expressed in terms of Bessel functions. This became the most widely-used model 
to describe the magnetic field structure of MCs. MCs are believed to represent 
about one third of the total number of ICMEs (Gosling, 1990; Richardson and 
Cane, 2010). However, a number of other studies proposed that a larger fraction 
of ICMEs are MCs; for example, the study by Li et al. (2011) and Marubashi 
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(2000) estimated the proportion of MCs to be 50% and 80%, respectively. Other 
authors have studied the variation of the proportion of MCs among ICMEs 
during the solar cycle (Kilpua et at, 2011; Richardson and Cane, 2010) and 
found that it varies from close to 100% near solar minimum to about 20-25% 
near solar maximum. 

2. Ejecta or irregular and weak magnetic field ICMEs: These configurations 
have some but not all the properties of MCs. Typically, their magnetic field 
is too weak and/or irregular or without much rotation. Some of these ejecta 
can be further classified into sub-categories, for example, "magnetic cloud-like" 
(Zhang et at, 2007) or magnetic flux rope (Rouillard, 2011). It has been recently 
proposed that some of these events are actually MCs whose legs passed through 
the observing spacecraft (Marubashi and Lepping, 2007), which has been con- 
firmed with STEREO observations (Mostl et al, 2010). According to this type 
of studies, ejecta are MCs but are not observed as such due to an observational 
bias. 

3. Complex ICMEs (Complex ejecta), which may result from the interaction 
of successive CMEs or from the interaction of CMEs with complex solar wind 
structures and streams (Burlaga, Plunkett, and St. Cyr, 2002; Wang, Ye, and 
Wang, 2003; Lugaz, Manchester, and Gombosi, 2005). Double rotation within 
one MC have also been reported from Ulysses (Rees and Forsyth, 2004) and 
ACE observations (Steed et al., 2011), showing that even isolated ICMEs may 
produce complex ejecta. 

Understanding the nature of the magnetic field in ICMEs (magnetic clouds 
and ejecta) is a crucial part of comprehending their behavior and the way they 
affect the interplanetary medium and Earth' s magnetosphere throughout their 
propagation. Thus there has been several attempts to describe the structure 
of the magnetic field in ICMEs using models and magnetic field reconstruction 
codes. Figure 1 gives some examples of proposed models. 

In order to explain the behavior of the magnetic field in MCs, early models 
were developed to accommodate Burlaga et al. (1981) 's definition of a MC. Most 
of these models restrict the internal structure of the magnetic field to the force- 
free configuration where the magnetic field, B, is described as V x B = aB. 
Several magnetic field reconstruction and fitting models have been built upon 
this definition; among the most used is the model by Burlaga (1988) and Lepping, 
Burlaga, and Jones (1990) which assumes a cylindrically symmetric solution and 
a constant a across the cloud. 

Marubashi (1986) also adopted a force-free model, but without assuming a 
constant a. Farrugia et al. (1993) and Farrugia, Osherovich, and Burlaga (1995) 
introduced a cylindrical model incorporating self-similarly expansion in two 
initially force-free constant-a codes, a cylindrical and a spherical models, and 
noted that the cylindrical configuration did not maintain the force-free state 
after it starts to expand. Thereafter, Shimazu and Vandas (2002) provided a 
modification of the mathematical formalism in order to keep the self-similarly 
expanding cylindrical model as force-free. 

Several other models have emerged thereafter which make different assump- 
tions: fitting to a cylinder of elliptical cross section by Hidalgo, Nieves-Chinchilla, 
and Cid (2002) and Hidalgo (2003); fitting to a kinematically distorted flux 
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Figure 1. An overview of different models to understand MCs. (a) Goldstein (1983) proposed 
that MCs are force-free structures, which was confirmed by Marubashi (1986) (figure adapted 
from Bothmer and Schwenn,1998) (b) Lepping et al. (1990) succeeded in fitting the Lundquist 
model to in situ magnetic field observations after an idea by Burlaga (1988). (c) Farrugia et 
al. (1993) included expansion into the model, (d) Schematic representation of the spacecraft 
trajectory inside of a MC of ellipitical cross-section following the model of Hidalgo et al., 
2002. (e) Hu and Sonnerup (2001) model magnetic flux ropes using magneto-hydrostatic 
reconstruction technique without a predefined geometry, (f) Marubashi and Lepping (2007) 
include curvature into the classic model. 



rope (Owens, Merkin, and Riley, 2006; Vandas and Romashets, 2003; Van- 
das et al., 2006; Dcmoulin and Dasso, 2009); non-cylindrical flux rope fitting 
(Mulligan and Russell, 2001; Owens et al., 2012); torus fitting (Romashets and 
Vandas, 2003; Marubashi and Lepping, 2007). The Grad-Shafranov reconstruc- 
tion technique (Hu and Sonnerup, 2002; Mostl et al., 2009) assumes a structure 
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in magneto-hydrostatic (MHS) equilibrium with an invariant direction, and uses 
the Grad-Shafranov equation to describe the magnetic field in the structure. 
Therefore, this model does not only use the magnetic field measurements but 
also measurements of the plasma pressure, and it is only applicable to a struc- 
ture possessing axial symmetry. The Grad-Shafranov technique has been tested 
using multi-spacecraft measurements (Liu et al., 2008; Kilpua et al., 2009; Mostl 
et al, 2009) and, recently, improvement to the algorithm have been presented 
by Isavnin, Kilpua, and Koskinen (2011). 

While most techniques reconstruct or fit ICMEs as some type of twisted flux 
ropes, this result is not necessarily well tested. It has recently been shown by 
Al-Haddad et al. (2011) that the Grad-Shafranov technique is designed such 
that it will always reconstruct a helical flux rope from a rotating magnetic field 
observed by a single-spacecraft, even if the magnetic field is not a helical flux 
rope in 3-D. This is precisely because Grad-Shafranov reconstruction assumes 
an invariance along the cloud axis (so-called 2.5 D method). As was shown in 
Al-Haddad et al. (2011), this is equivalent to assuming a helical solution from 
the beginning. The assumption of invariance along the cloud axis is also made 
in all the other fitting techniques used in this article. Other recent observations 
pointing towards a more complicated geometry of ICMEs at 1 AU have been 
made by Kahler, Krucker, and Szabo (2011) and Farrugia et al. (2011). 

Using 2-D and 3-D numerical simulations, it is possible to show that a CME 
initiated at the solar surface (for example using a twisted flux rope) evolves 
through its interaction with the solar wind into a typical ICME (Riley, Gosling, 
and Pizzo, 1997; Riley et al, 2003; Odstrcil and Pizzo, 1999; Manchester et at, 
2004; Chane et al., 2006; Shen et al., 2007). This idea has now been confirmed 
by STEREO observations (see, for example, Savani et al., 2010). In these sim- 
ulations, synthetic satellite measurements made at 1 AU show a typical ICME 
structure with a fast shock, preceding a dense sheath and an ejecta. Synthetic 
coronagraphic and hcliospheric images are also able to reproduce typical views 
of CMEs (Lugaz, Manchester, and Gombosi, 2005; Manchester et al., 2008; Riley 
et al, 2008; Lugaz et al, 2009; Odstrcil and Pizzo, 2009). 

Riley et al. (2004) made a comparison of magnetic field reconstruction and 
fitting models for ICMEs, by fitting two different time series of a simulated 
ICME to five different techniques: three force-free models, the elliptical model 
by Hidalgo, Nieves-Chinchilla, and Cid (2002) and the Grad-Shafravov recon- 
struction code. The results of the reconstructions were then compared back to the 
2.5-D MHD simulation. One of the difficulties that the authors encountered was 
choosing the boundaries. Another one was recognizing a MC within the provided 
time series. Specifically, when the minimum variance technique is used and the 
impact parameter is large, it is hard to observe the properties characterizing a 
MC. 

Three of the models used in the study by Riley et al. (2004) are also used in 
the present study. We also added the self-similarly expanding cylindrical fitting 
of Marubashi and Lepping (2007) and only kept one classical -no expansion- 
force-free model. However, it is important to stress that in this project we apply 
the comparison to a list of observed events. Here, the authors address some of 
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the issues discussed in Riley et al. (2004) by performing more detailed analyses 
and proposing fixed boundaries. 

In this study, we compare the reconstruction of four different magnetic field 
fitting techniques for most of the fast ICMEs and MCs observed during solar 
cycle 23, for which a source region can be identified (see details in section 2). 
Our goals are to acquire a broader understanding of the magnetic field structure 
in MCs and to investigate if there are statistical differences between the codes 
for MCs and non-MC ejecta, and take steps towards deciphering if all ICMEs 
can be expressed as flux ropes if a sophisticated enough model is used. We do 
this by comparing the results of the different fitting techniques. In section 2, we 
describe the data used in this study followed by a succinct description of the 
different codes. In section 3, we discuss the full results for the 59 selected events 
and in section 4, we present two events for which we used the same boundaries 
for all the codes. We conclude and discuss our results for the nature of magnetic 
fields in ICMEs in section 5. 



2. Data and Models Overview 



2.1. CDAWData 



The data is taken from the list of Coordinated Data-analysis Workshop (CDAW) 
of MCs and ejecta during solar cycle 23. The list can be found at http:// 
cdaw.gsfc.nasa.gov/meetings/2010_fluxrope/LWS_CDAW2010JCMEtbl.html. This 
list contains all shock-driving ICMEs from solar cycle 23 for which it was possible 
to determine the source region, and for which the source region was within ±15° 
of disk center. This list is a subset of the list of all interplanetary shocks observed 
during solar cycle 23 from Gopalswamy et al. (2010). The list contains 59 ICMEs 
classified as follows: 24 events as Magnetic Clouds (MCs) and 35 events as Ejecta 
(EJ). 



2.2. Grad-Shafranov Reconstruction Technique 



For the Grad-Shafranov (GS) reconstruction, we use the code by Hu and Son- 
ncrup (2002) and followed the guidelines for using it discussed by Mostl et al. 
(2009). This is based on magneto- hydrostatic equilibrium of a system with an 
invariant direction and is a solution for what is basically a numerical boundary 
problem. The GS equation is: 



8 2 A 3 2 A _ dPt(A) 

where A(x,y) is the magnetic vector potential, P t the transverse pressure 
defined as P t = p + B 2 /(2/j,q) with p being the thermal pressure. The total 
magnetic field is given by B = VA x e z + B z e z . The first step to solve this 
equation is to determine the invariant axis, z. This is done first using a minimum 
variance analysis on the magnetic field components and by finding a frame where 



SOLA: Comparison_PaperAfter_revs.tex; 1 October 2012; 0:39; p. 6 



CDAW-II: CMEs Magnetic Field Fitting and Reconstruction Codes 



the transverse pressure is a single- value function of the magnetic vector potential, 
A. After this, the GS equation can be solved numerically. 

The major assumptions of the technique are: the structure is assumed to 
be invariant along the ICME axis (so called 2.5 D), and it is assumed to be 
time- independent during the whole measurement by the spacecraft. Compared 
to other techniques, it does not include expansion or a toroidal geometry. One 
of the potential issues with the Grad-Shafranov technique is the stability of 
the integration and the influence of the solver and stabilization procedure on 
the results. However, Grad-Shafranov reconstruction has been tested by using 
observations of ICMEs by two spacecraft at a significant separation with respect 
to the ICME size (Liu et al., 2008; Mostl et al, 2009; Mostl et al, 2009). Some 
of the stability issues have also been recently addressed in Isavnin, Kilpua, and 
Koskinen (2011). 

One of the main advantages of the GS method is that the shape of the flux 
rope's cross-section and the number of flux ropes inside the reconstructed interval 
are both not pre-defined, but are an output of the technique. Another advan- 
tage is that the boundaries of the ICMEs are chosen through an optimization 
procedure to make the P t (A) function as close to single- valued as possible. The 
Pt(A) function is fitted with a polynomial function with exponential tail(s) and 
only these events for which the fitting residue is small (typically less than 0.1) 
and the fitting appears visually correct arc deemed as successfully reconstructed. 
Further details can be found in Hu et al. (2004). 

GS reconstruction gives the following results: flux rope orientation (longitude 
and latitude), impact parameter (closest approach of the spacecraft to the MC 
axis), magnetic flux (axial and poloidal), axial current, the cross-section shape 
and the central field strength. The GS reconstruction feasibility was checked for 
every event in the CDAW list. We used primarily data from the Wind spacecraft 
(MFI (Lepping et al, 1995) and SWE (Ogilvie et al, 1995) instruments), and 
switched to ACE (MAG (Smith et al, 1998) and SWEPAM (McComas et al, 
1998) instruments) if there was a data gap at Wind. The GS reconstruction 
method was successful for: 

• full list: 20 events out of 59 (34 %) 

• 17 MC events out of 24 (71 %) 

• 3 EJ events out of 35 (8 %) 

This means that the applicability of using GS reconstruction is closely tied to 
the definition of a magnetic cloud (higher-than-average field strength, smoothly 
rotating magnetic field vector, low j3, low T p ). Typically, this means that, for 
these cases, one magnetic field component (B y or B z ) needs to be bipolar and 
the other unipolar, depending on the orientation of the flux rope (see discussion 
in Bothmer and Schwenn, 1998). Also, we were able to reconstruct about 1/3 of 
the events in the complete list, which is similar to the often quoted ratio between 
all ICMEs and those which contain parts that satisfy the MC definition. 

2.3. Sclf-similarly Expanding Cylinder Fitting Technique 

This model is a force-free model that assumes a spatially constant a and performs 
the fitting in cylindrical coordinates. Here, a is chosen fixed for every event 
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Figure 2. Axial magnetic field for the Grad-Shafranov (red diamonds) and self-similarly 
expanding cylinder (purple stars) models as compared with the axial magnetic field for the 
FFCA model. The red and purple dashed lines show the value of the best-linear fit to the 
Grad-Shafranov and self-similarly expanding cylinder models, respectively. 



and equal to 2.4. The model takes into account self-similarly expansion, as was 
initially suggested by Farrugia et at, Farrugia et al. (1992, 1993). The original 
inclusion of the self-similarly expansion has been found to cause a deviation from 
the force-free state during expansion. Some refinements proposed by Shimazu 
and Vandas (2002) have been incorporated into the model to preserve the force- 
free state for the duration of the MC travel. As a result, the magnetic field 
intensity B and the radius r of the fitted MC at a given time as produced by 
the model are defined as follows: 



-Bo 



B= „ ■ 1,2 ( 2 ) 



r = r (l + ~) (3) 

^exp 

where Bq is the magnetic field intensity on the axis of the CME cylinder, t cxp 
is the expansion time, t is the time since the first encounter of the spacecraft 
with the CME, ro is the distance between the cylinder axis and the spacecraft 
at the time of the first encounter. 

The model also provides values for the longitude and latitude of the axis, as 
well as the impact parameter, which is defined to as the distance between the 
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Figure 3. Same as Figure 2 but for the longitude and latitude of the axis of the reconstructed 
ICME. The yellow squares are for the elliptical cross-section fitting model. The dashed lines 
highlight to angles ±45° from the FFCA angle. 

closest axis of cylinder and the path of the spacecraft. The model fits simul- 
taneously the following quantities for both chiralities (positive and negative): 
Bo, r , t cxp , the longitude, latitude and impact parameter. The final selection is 
made by choosing the sign of the chirality, which gives the best fit. The details 
of the model can be found in Appendix A of Marubashi and Lepping (2007). 
The self-similarly expanding cylinder fitting method was successful for: 

• full list: 50 events out of 59 (85%) 
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• 22 MC events out of 24 (91.5 %) 

• 28 EJ events out of 35 (80 %) 

There were seven events, which were found to be contracting (negative i cxp ) 
and three events, for which the expansion time was greater than 500 hours, which 
correspond to less than 10% change in the axial magnetic field strength over the 
crossing of a typical magnetic cloud (<~ 1 day). 

2.4. Force-Free Constant-Alpha Reconstruction 

This force- free fitting model is based on the force- free flux constant-alpha (FFCA)| 
rope developed by Burlaga (1988) and then optimized by Lepping, Burlaga, and 
Jones (1990). In our model, a is fixed to equal 2.4. This implies the flux rope's 
magnetic field on the outer shell is almost completely poloidal. The value of 2.4 
is not optimized and it could be changed (i.e. reduced in value) if it is believed 
that on the outer shell of the flux rope the field lines are still helical. This might 
be the case if the rear or front edge of the flux has been "eroded" away from 
its ideal case (perhaps by reconnection) as proposed by Dasso et al. (2006) and 
recently confirmed by Ruffenach et al. (2012). The a parameter has also been 
used as a free variable is some models (e.g. Lepping, Burlaga, and Jones, 1990). 
Since a is constant and fixed to equal 2.4 for the FFCA and for the self-similarly 
expanding cylinder models, it is clear that these two models are identical when 
the expansion is equals to 0. 

The flux rope axis is calculated by Minimum Variance Analysis (MVA) . MVA 
is a technique originally developed for solving the normals to a tangential discon- 
tinuity (TD) (Sonnerup and Scheible, 1998). In order to obtain the orientation 
of an ICME, the axis is determined by the intermediate eigenvector, as pointed 
out by Goldstein (1983). 

The chirality of the flux rope is determined by the sense of direction the 
magnetic field rotates in the MVA reference frame. The impact parameter and 
magnetic field magnitude are the only two free variables. They are optimized 
using a non- linear optimization routine (Nelder and Mead, 1965). For each event, 
the fitting was attempted with two different boundaries: those from Lepping 
et al. (2006) and those selected by K. Marubashi. In addition, for eight events, 
a different set of boundaries were chosen directly by N. P. Savani as the other 
two sets of boundaries were judged as not satisfactory for these events. Overall, 
the method was successful for: 

• full list: 55 events out of 59 (93%) 

• 24 MC events out of 24 (100%) 

• 31 EJ events out of 35 (88.5%) 

2.5. Elliptical Cross-section Fitting Model 

The elliptical, non-force-free model (Hidalgo, Nieves-Chinchilla, and Cid, 2002; 
Hidalgo, 2003) considers that ICMEs have elliptical cross-section due to their 
interaction with the solar wind. In this model, the expansion of the local cross- 
section is taken into account. Two characteristics distinguish this model from the 
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others used in this paper: 1) non-force-free condition is assumed, and, 2) an ellip- 
tical cylindrical coordinate system is chosen to solve the Maxwell equations. The 
model assumes a flux rope magnetic field geometry and treats the radial and axial 
electric currents separately. The local solutions of the Maxwell equations provide 
the three magnetic field components to be transformed into spacecraft coordi- 
nates. The only assumptions are that the radial and axial components of the 
current density are constant. The cross-section of the flux ropes is not required 
to be circular. This model was compared with the circular cross-section model of 
Lepping, Burlaga, and Jones (1990) in Nieves-Chinchilla, Hidalgo and Sequeiros 
(2005) and compared with information from remote-sensing measurements in 
Nieves-Chinchillaef al. (2012). 

The ellipticity, the current density components, expansion time of the cloud 
and the axial magnetic field are determined by fitting the magnetic field data as 
explained in section 2 of Hidalgo (2003) . Multiple regression analysis is used to 
infer the spacecraft trajectory through the flux rope. The method converged for 
all events. However, the fit goodness was reported with a flag: 0=bad; l=good; 
and, 2=very good. Such "goodness" is based on the correlation coefficient, fitting 
residue and visual goodness. In the rest of the analysis, we take into account 
just fitting results with quality 1 or 2 (i.e., removing bad fits). In that case, the 
elliptical cross-section fitting method was successful for: 

• full list: 39 events out of 59 (66 %) 

• 24 MC events out of 24 (100 %) 

• 15 EJ+ events out of 35 ( 43%) 

2.6. Selection of the ICME's Boundaries 

Selection of the boundaries of an ICME is still a challenging problem. This 
is because there is no single definitive characteristic of ICMEs. Some authors 
base the ICME boundaries purely on magnetic properties; others use low proton 
temperature as one of the criteria; while others use charge state information. 
Many of the in situ signatures have been taken into account in order to determine 
the boundaries of ICMEs. However, each of those indicators provide different 
boundaries which leaves the problem still unsolved. In this study, we let each 
group determine the boundaries for the ICMEs. For three events, we re-run 
the analysis using pre-set boundaries as given by the self-similarly expanding 
cylinder fitting technique. 

3. Comparative analysis of the fitting parameters Results 

Out of 59 events on the CD AW list, the reconstruction of the magnetic field 
was successfully done for 20 events, mostly MCs, by the GS technique. By 
comparison, the FFCA fitting and the self-similarly expanding cylinder methods 
were able to fit more than 85% of the list, and the elliptical cylinder fitting 
method 66%. In the following sections, we compare the fitted parameters from 
different codes with that from the FFCA code because (i) it is the simplest 
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code in our study, and (ii) it succeeded in fitting most of the cases with the 
least deviation from the results of the other models. The number of cases that 
have been analyzed in common between the FFCA code and the other codes is 
as follows: 20 with the GS reconstruction technique, 45 with the self-similarly 
expanding cylinder fitting method, and 37 with the elliptical cylinder fitting 
method. 

3.1. Chirality 

For all 20 events reconstructed by the GS code, the chirality (sign of the magnetic 
helicity) was the same as for the self-similarly expanding cylinder and the FFCA 
model. There are three events for which the FFCA and self-similarly expanding 
cylinder found opposite signs for the helicity (13, 15 and 52). This shows that 
even the simplest result of fitting and reconstruction methods is not necessarily 
well constrained. The difference is primarily due to the choice of boundaries. In 
two of these three cases (13 and 15), when the FFCA fit is performed again with 
a different choice for the boundaries, it returns the same sign of helicity as that 
found by the self-similarly expanding cylinder. In fact, it is worth noting that 
the chirality found with the FFCA differs for 7 events depending on whether the 
boundaries are chosen following Lepping et al. (2006) or that chosen specifically 
for this study by K. Marubashi. 

One event is more ambiguous. For event 52, the FFCA fit is relatively poor 
and gives a left-handed ejecta for the boundaries selected by K. Marubashi as 
well as for those from Lepping et al. (2006) , whereas the self-similarly expanding 
cylinder fit gives a right-handed ejecta. It should be noted that the best-fit axial 
magnetic field by FFCA for this event is less than 6 nT for the two sets of 
boundaries, reflecting a very weak ejecta. 

This overall result about the chirality may appear inconsistent with what 
was found in previous studies by Dasso et al. (2003) and Dasso et al. (2006) 
where the magnetic flux and helicity was found to be nearly independent of the 
fitting and reconstruction model, for two well observed MCs. Here, we find that 
it is generally true, consistent with these results. However, we also show that 
for ejecta, for which the boundaries are not well defined because of incomplete 
rotation, fittings with different boundaries can yield opposite chirality for the 
same event. 

3.2. Maximum Magnetic Field Strength 

For all the common cases between the FFCA code and one of the other codes in 
this study, we calculated the correlation of the fitted axial magnetic fields. For the 
self-similarly expanding cylinder fitting method, we use the axial magnetic field 
halfway through the cloud, corresponding approximately to the magnetic field 
at the closest approach to the cylinder's axis, as calculated from Equation (2). 
We do not compare the value of the axial magnetic field for the elliptical cross- 
section fitting model with that of the other models. This is because the axial 
magnetic field in this model is the sum of two fields: the magnetic field created 
by the currents and the axial field which is a parameter from the model. Figure 2 
summarizes the magnetic field results. 
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Overall, the correlation between the FFCA, GS and self-similarly expanding 
cylinder code is quite good (~0.84). The axial magnetic field value from the 
GS code tends to be slightly higher than that from the FFCA especially for 
stronger magnetic fields, whereas the self-similarly expanding cylinder code tends 
to return higher values for the magnetic field for weak magnetic fields. 

3.3. Orientation of the ICME Axis: Longitude and Latitude 

The orientation of the ICME as it arrives at Earth is an essential parameter, since 
it is often used to determine the CME rotation as it propagates in the heliosphere 
(Yurchyshyn et ai, 2001). In addition, knowing precisely the orientation of an 
ICME is particularly important since the orientation of an ICME is known to 
be related to its geo-effectiveness (see, for example Zhao and Hoeksema, 1998). 
As can be seen in Figure 3, the orientation of some ICME events differs greatly 
depending on which model is used. In fact, only for one CME (event #44) did 
all events give an orientation of the axis within ±45° of each other. 

First, we focus on a comparison of the orientation between the Grad-Shafranov,| 
self-similarly expanding cylinder and FFCA models for the 20 events successfully 
reproduced by all three methods. For seven events, the three methods give 
an orientation within ±45° of each other. For eleven other events, the Grad- 
Shafranov reconstruction technique gives an orientation of the axis within ±45° 
of one of the other two methods (but not the other). In fact, although the FFCA 
and self-similarly expanding cylinder fitting methods are closely related, the 
axis of the fitted cylinder is found to be in agreement for only about 45% of the 
events (21 out of the 45 common cases; for 14 of these, the orientation is within 
±30°). While this may appear surprising considering the similarities between the 
two methods, it should be noted that the force-free method without expansion 
relies on the minimum-variance analysis to obtain the orientation of the flux 
rope, whereas the expanding cylinder method relies on a multi-parameter fitting 
procedure to minimize the deviation between the model and the measurements. 
In addition, for many events, both methods used different boundaries. 

Comparing the direction obtained from the elliptical cross-section model with 
that from the force- free cylindrical models (with and without expansion), we 
find that the direction from all three models only agree for one event. For seven 
additional events, the direction from the elliptical cross-section model is within 
±45° of the direction from the FFCA fitting model; for three additional cases, 
the same is true as compared to the direction from the self-similarly expanding 
cylinder model, and for one event it is true as compared with the direction for 
the Grad-Shafranov reconstruction technique. 

Overall, there are only 34 events reconstructed by three methods or more. 
For 30 of these events, the orientation of the axis of the reconstructed ICME 
is consistent between at least 2 methods. The same is only true for three out 
of the 18 events reconstructed by two methods only. Finally, seven events were 
reconstructed by one method only. 

3.4. Difference Between MCs and non-MC Ejecta 

Here, we briefly compare the average values of some parameters for the MCs 
with those for the non-MC ejecta. We only focus on the FFCA and self-similarly 
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expanding cylinder, because the Grad-Shafranov reconstruction methods is only 
able to reconstruct three non-MC ejecta (out of 35) and the elliptical cross- 
section fitting has different set of parameters than these two models. Consistent 
with their definitions, non-MC ejecta are found to have, on average, weaker 
magnetic field than MCs: 13.75 nT vs. 22.1 nT (FFCA) and 17.7 nT vs. 34.4 nT 
(expanding cylinder). We also find that the average impact parameter is larger 
for non-MC ejecta than for MCs: 0.73 vs 0.54 (expanding cylinder) and 0.32 vs. 
0.29 (FFCA). This may indicate that ejecta are indeed flux rope-like magnetic 
clouds but they are measured at a large impact parameter, resulting in weaker 
magnetic field and less well-defined properties. Note that the average impact 
parameter for the events reconstructed by the GS method is 0.29 (FFCA) and 
0.47 (expanding cylinder), once again showing that the GS method was only 
successful for MCs or MC-like ejecta. 

4. Detailed Analyses of Three Events 

We focus on events 10, 20 and 27 which we fitted using the same boundaries. 
These were those chosen for the self-similarly expanding cylinder fitting method. 

4.1. Event 10 

This event is a left-handed non-MC ejecta observed in 1999 June 26-28. This was 
one of the events resulting in a "poor" fit for the elliptical cross-section model and 
it was not initially successfully reconstructed by the Grad-Shafranov method. For 
this method, one problem was the presence of a shock propagating inside the 
ICME at 19:16 UT on June 26. Because the GS reconstruction method attempts 
to find the best orientation of the ICME by requiring the function Pt(A) to be 
single- valued, the presence of a shock in the magnetic and plasma pressure in one 
side of the cloud makes this an impossible task. This has been previously noted 
in Kilpua et al. (2009) and Mostl et al. (2009). It is not a problem encountered 
by the other fitting methods. 

Initially, some researchers found two possible flux rope intervals: one shortly 
after the shock wave on June 26 lasting about 16 hours (referred hereafter as 
10a) and one corresponding to the initial interval fitted with the self-similarly 
expanding cylindrical code of 6.5 hours at the end of June 27 and in the beginning 
of June 28 (referred hereafter as 10b from 20:30 UT on June 27 to 03:00 UT on 
June 28). We only discuss here the results from the longer, first interval. The 
intervals found by different researchers are shown in Table 1. 

The analysis was performed again for the interval 10a using for interval June 
26 from 06:00 to 19:00. In Figure 4, we show the fitting results for the differ- 
ent methods corresponding to this interval. Using the same interval for all the 
methods, the magnitude of the magnetic field in the ICME is found to match 
relatively well between the two force-free and the Grad-Shafranov methods and 
to correspond to a value ~ 15 — 20 nT. The orientation of the cloud is found to be 
consistent with a latitude of 0° ± 12° and a longitude of 65° ± 25° between these 
three methods, whereas the elliptical cross-section model finds a value about 90° 
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away. This example illustrates the importance of the choice of boundaries for the 
study of magnetic clouds and ICMEs, as three of the four codes find very similar 
parameters when the same boundaries are chosen. 

Table 1. Intervals used by different researchers for the events 10, 20 and 27. In bold are the 
intervals chosen when the events were compared with the same boundaries. In parentheses 
are the day of the month corresponding to the interval boundary. Event 10, 20 and 27 took 
place in June 1999, July 2000 and November 2000, respectively. 



Event 


10 (June) 


20 (July) 


27 (Nov.) 


GS 


N/A 


08:36(27)-19:24(27) 


23:56(06)-17:20(07) 


FFCA 


06:00(26)-22:00(26) 


07:30-19:30 


22:10(06)-17:10(07) 


Exp. Cylind. 


06:00(26)-19:00(26) 


07:30-19:30 


22:10-17:10 


Ellip 


22:36(26)-02:36(27) 


10:10-20:40 


23:08-17:45 



4.2. Event 20 

This event is a left-handed non-cloud ejecta observed in 2000 July 27. It is 
one of the three non magnetic cloud ejecta successfully reconstructed using the 
Grad-Shafranov reconstruction method. Initially, as we let each researcher select 
the ejecta's boundaries, the direction of the cloud axis was in good agreement 
between the four methods: a longitude of 240° ± 30° and a latitude of 22° ± 12°. 
All methods were within 45° of each other, except between the self-similarly 
expanding cylinder and the elliptical cross-section cylinder which were about 
51° apart. 

The impact parameter was found to be close to for the Grad-Shafranov but 
larger than 0.5 for the other three methods. The magnitude of the axial magnetic 
field was about 20 nT for the Grad-Shafranov method and about 8 nT for the 
two force-free constant-alpha models. 

As can be seen from Table 1, all the different researchers chose boundaries 
very close to each other and choosing exactly the same boundaries does not 
drastically change the results, except that the maximum axial magnetic field for 
the Grad-Shafranov model was 12 nT with an impact parameter of about 0.25, 
more consistent with the values found by the other models. In Figure 5, we show 
the results of the fit using the boundaries chosen for the self-similarly expanding 
cylinder fit. 

4.3. Event 27 

This event is a left-handed magnetic cloud observed in 2000 November 6-7. All 
methods agree on the maximum magnetic field strength (<~ 25 nT for Grad- 
Shafranov and the FFCA model and ~31.5 nT for the self-similarly expanding 
cylinder). However, there is initially little agreement regarding the orientation of 
the magnetic cloud with a latitude of 14° ± 14° and a longitude of 145° ± 55° for 
the four methods. As for the previous event, the boundaries chosen independently 
by different researchers were relatively consistent. One significant difference is the 
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inclusion or not of a region of strong magnetic field without much rotation at the 
"back" of the MC (until 02:40 on November 8). Previous studies (Dasso et al., 
2006; Dasso et al, 2007; Mostl et al, 2008) have discussed how reconnection 
between a MC and the ambient solar wind during the ICME propagation may 
result in this type of regions which, actually, belong to the MC. Each researcher 
decided exactly where to end the MC. Since there is no counterpart flux on the 
MC front to the flux on the back of the MC, it is probably better not to take 
this region into account. 

Using the same boundaries, there is a better agreement that the magnetic 
cloud has a low inclination (latitude ~ 5) but not much better agreement 
about the precise orientation or the value of the impact parameter. The fits 
and reconstructions are shown in Figure 6. 

5. Summary and Conclusions 

In this article, we have presented a comparison of the reconstruction results of 59 
ICMEs measured in situ during solar cycle 23 using four different reconstruction 
or fitting methods. The events our study focused on were chosen because their 
source region was known and was within ±15° of the central meridian. Our data 
included 24 events identified as magnetic clouds (MCs) following the definition 
of Burlaga et al. (1981) and 35 non-MC ejecta. We find that the two force-free 
techniques (with or without expansion) are able to fit the vice majority of the 
events (>85 %) including more than 75% of the ejecta. The Grad-Shafranov re- 
construction method which assumes magneto-hydrostatic equilibrium and makes 
use of the plasma pressure is found to only work for clear magnetic clouds, as 
it was only successful in reconstructing 8% of the non-cloud ejecta (3/35). In 
addition, the Grad-Shafranov reconstruction is not able to reconstruct ICMEs 
into which a shock is propagating as it violates its assumption. However, one of 
the advantages of this method is the fact that the boundaries of the ICME are a 
result of the reconstruction, and they do not need to be chosen before performing 
the fits, as for the other methods. The expanding elliptical cross-section was able 
to fit with a reasonable result about two-thirds of the events (all magnetic clouds 
and about 40% of the ejecta). 

We have found that even the ICME chirality (sign of the magnetic helicity), 
one of the simplest possible reconstruction parameters, is not necessarily well 
constrained and opposite chirality can be found for the same event because of 
differences in the choice of the boundaries. In most cases, the choice of boundaries 
alone explain the difference between different methods. 

Regarding the maximum magnetic field strength, we have found that, typi- 
cally, the results between Grad-Sharanov method, FFCA and self-similarly ex- 
panding fittings are well correlated and within about 25%. The Grad-Shafranov 
method returns larger values than the FFCA fitting method for strong magnetic 
fields, whereas the force-free model which includes self-similar expansion returns 
higher values than the FFCA fitting method for weak magnetic fields. Gulisano 
et al. (2005) previously compared the maximum magnetic field at the cloud axis 
for 20 magnetic clouds using four fitting methods and found that the variability 
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between methods for the same event was significantly less than the variability 
between events for the same method (see their figure 4). 

For only one event, all four methods find an orientation of the ICME axis 
within ±45°. Considering only the two force- free fitting methods (with and 
without expansion) and the Grad-Shafranov reconstruction method, this is the 
case for about 30% of the ICMEs (7 out the 20 common cases). Directly com- 
paring the two force- free constant-alpha codes (with and without expansion), 
they return ICME axis within ±30° in only <~ 30% of the cases and within 
±45° in only ~ 45% of the cases. While these two methods assume very similar 
geometries for the ICME, the fitting procedures are very different: minimum 
variance analysis is used to determine the axis' direction for the FFCA code, 
whereas a fitting procedure is used for the expanding cylinder method. This is 
likely to account for the difference in the direction of the ICMEs' axis. Finally, 
the elliptical cross-section model typically gives different orientation as compared 
to the other methods, being in agreement with at least one other method for 
only about 30% of the events. These results were found by letting researchers 
select their own boundaries. 

It should be noted, that when only two methods successfully reconstruct an 
event, the orientation found by these two methods usually disagree (by more 
than ±45° in 15 out the 18 such cases). When the two force- free constant- 
is fitting methods and the Grad-Shafranov reconstruction method successfully 
reconstruct an event, the orientation from the Grad-Shafranov method is within 
±30° of that from one of the force-free constant-a method in 70% of the cases 
(14 out of 20 such events). Combining all four methods, more than 65% (23/34) 
of the events reconstructed by three or more methods have at least two methods 
giving orientation within ±30° . It is therefore clear that having multiple methods 
able to successfully fit or reconstruct the same event gives more reliable results 
regarding the orientation of the ICME axis. 

We further quantified the importance of the selection of the boundaries by 
performing the analyses again for three events using the same boundaries for all 
methods. By using the same boundaries, we find a better agreement between the 
different codes regarding the ICME magnetic field strength, and, for one event, 
a better agreement for the ICME orientation. 

Finally, for two methods, we have compared the impact parameter and mag- 
netic field strength of MCs and non-MC ejecta. We find that the non-MC ejecta 
have, as expected, weaker magnetic field but we also find some evidence that 
the impact parameter is, on average, larger for non-MC ejecta than for MCs. 
This is a first statistical hint that most (or all) ICMEs observed at 1 AU may in 
fact, be magnetic flux ropes, but depending on how they impact the observing 
spacecraft, they may or may not be recognized as magnetic clouds. Having 
performed fitting and reconstructions for 59 ICMEs with four different codes, 
we will further investigate the nature of the magnetic fields in ICMEs as well as 
how to further improve reconstruction methods in the future. 
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CDAW-II: CMEs Magnetic Field Fitting and Reconstruction Codes 




Figure 4. Fitting for event 10a by the different codes. Top left: Self-similar expanding cylinder 
code fit shown in red. The different panels are the total magnetic field, B x , B y , B z , the magnetic 
field fluctuation, the solar wind speed, proton density, alpha-to-proton ration, temperature and 
plasma f3 from top to bottom. Bottom left: FFCA code fit shown in purple. The panels show 
from top to bottom, the radial, tangential, normal (R, T and N) components, the longitude 
and latitude angles (in RTN) and magnitude of the magnetic field, the proton velocity, density 
and temperature. Top right: GS reconstruction map, the axial magnetic field is color-coded 
and the black contours are magnetic field lines in the plane of the cross-section of the 1CME. 
Observed magnetic field components in this plane are shown with white vectors and velocity 
in green. Bottom left: Elliptical cross-section fit in red. The different panels show the magnetic 
field, its three components, the plasma j3 and the solar wind speed from top to bottom. 
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Figure 5. Same as Figure 4 but for event 20. 
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Figure 6. Same as Figure 4 but for event 27 
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